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- - - ^ We performed Self-Consistent Greens Function (SCGF) calculations for symmetric nuclear matter 

^y , using realistic nucleon-nucleon (NN) interactions and effective low-momentum interactions [Viow-k), 

>-^ ' which are derived from such realistic NN interactions. We compare the spectral distributions re- 

sulting from such calculations. We also introduce a density-dependent effective low-momentum 
interaction which accounts for the dispersive effects in the single-particle propagator in the medium. 
U 
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I. INTRODUCTION 



K*" ' The description of bulk properties of nuclear systems starting from realistic nucleon-nucleon (NN) interactions is a 

^^ • long-standing and unsolved problem. Various models for the NN interaction have been developed, which describe the 
pZ. \ experimental NN phase shifts up to the threshold for pion production with high accuracy[2l3,l2, ^4J- A general feature 
_^ . of all these interaction models are strong short-range and tensor components, which lead to corresponding correlations 
f— ^ ' in the nuclear many-body wave-function. Hartree-Fock mean-field theory, which represents the the lowest-order many- 
\^ , body calculations one can perform with such realistic NN interactions, fails to produce bound nuclei 0, Q precisely 
^^ ■ because Hartree-Fock does not fully incorporate many-body correlation effects. 

That correlations beyond the mean field are important is supported by experiments exploring the spectral distribu- 
tion of the single-particle strength. One experimental fact found in all nuclei is the global depletion of the Fermi sea. 
A recent experiment from NIKHEF puts this depletion of the proton Fermi sea in -^^^Pb at a little less than 20% [3 
y . in accordance with earlier nuclear matter calculations [S|. Another consequence of the presence of short-range and 
tensor correlations is the appearance of high-momentum components in the ground state wave-function to compensate 
for the depleted strength of the mean field. Recent JLab experiments Q indicate that the amount and location of 
this strength is consistent with earlier predictions for finite nuclei [lOj and calculations of infinite matter [ll| . 

These data and their analysis, however, are not sufficient to allow for a detailed comparison with the predictions 
derived from the various interaction models at high momenta. In this paper, we want to investigate a possibility to 
Ci I separate the predictions for correlations at low and medium momenta, which are constrained by the NN scattering 
matrix below pion threshold, from the high momentum components, which may strongly depend on the underlying 
model for the NN interaction. For that purpose we will perform nuclear many-body calculations within a model space 
that allows for the explicit evaluation of low-momentum correlations. The effective Hamiltonian for this model space 
will be constructed from a realistic interaction to account for for correlations outside the model space. 

This concept of a model space and effective operators appropriately renormalized for this model space has a long 
history in approaches to the nuclear many-body physics. As an example we mention the effort to evaluate effective 
operators to be used in Hamiltonian diagonalization calculations of finite nuclei. For a review on this topic see e.g. |l2J . 
The concept of a model space for the study of infinite nuclear matter was used e.g. by Kuo et al. p^ll4Lll5l |. Also 
the Brucckncr-Hartree-Fock (BHF) approximation can be considered as a model space approach. In this case one 
restricts the model space to just one Slater-determinant and determines the effective interaction through a calculation 
of the G-matrix, the solution of the Bethe-Goldstone equation. 
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The effective hamiltonians for such model space calculations have frequently been evaluated within the Rayleigh- 
Schrodinger perturbation theory, leading to a non-hermitian and energy-dependent result. The energy-dependence 
can be removed by considering the so-called folded-diagrams as has been discussed e.g. by Brandow[i^ and Kuo[l7l|. 
We note that the folded-diagram expansion yields effective interaction terms between three and more particle, even 
if one considers a realistic interaction with two-body terms onlv[l3.ll9l|. 

During the last years the folded-diagram technique has been applied to derive an effective low-momentum potential 
Viouj-fcHfl from a realistic NN interaction. By construction, Viow-k potentials reproduce the deuteron binding-energy, 
the low-energy phase shifts and the half-on-shell T matrix calculated from the underlying realistic NN interaction 
up to the chosen cut-off parameter. The resulting Vioiu-k turns out to be rather independent on the original NN 
interaction if this cut-off parameter for the relative momenta is below the value of the pion-production threshold in 
NN scattering. The off-shell characteristics of the Viow-k effective interaction are not constrained by experimental 
data and can influence the many-body character of the interaction. 

For finite nuclei we find that one does indeed obtain different binding energies for ^^O depending on the underlying 
NN interaction from which one derives the Viow-k interaction. For example, using coupled-cluster techniques at the 
singles and doubles level (CCSD) |21j we find binding energies for ^^O at a lab-momentum cutoff of A = 2.0 fm"-'^ to 
be -143.4 ± 0.4 MeV and -153.3 ± 0.4 MeV for the N^LO j^] and CD-Bonn two-body interactions, respectively The 
CCSD calculations were carried out at up to 7 major oscillator shells (with extrapolations to an infinite model space) 
using the intrinsic Hamiltonian defined as H = T — Tcm + Viow-k where Tcm is the center of mass kinetic energy. 

Attractive e nerg ies arc obtained if such a Viow-k interaction is used in a Hartree-Fock calculation of nuclear matter 
or finite nuclei|22l l23l|. High-momentum correlations, which are required to obtain bound nuclear systems from a 
realistic NN interaction (see above) are taken into account in the renormalization procedure which leads to Viow-k- 
Supplementing these Hartree-Fock calculations with corrections up to third order in the Goldstone perturbation theory 
leads to results for the ground-state properties of ^^O and ^^Ca^ which are in fair agreement with the empirical data|23|. 
(One should note that Tcm was not included in these calculations.) Calculations in infinite matter demonstrate that 
Viow-k seems to be quite a good approximation for the evaluation of low-energy spectroscopic data. The results for 
the pairing derived from the bare interaction are reproduced |23l| . T he p rediction of pairing properties also agree with 
results obtained phenomenological interactions like the Gogny force[2J,|23. The Viow-k interaction also yields a good 
approximation for the calculated binding energy of nuclear matter at low densities. 

At high densities, however, BHF calculations using Viow-k yield too much binding energy and do not reproduce 
the saturation feature of nuclear matter |23i]. This is due to the fact that Vjom-fc does not account for the effects of 
the dispersive quenching of the two-particle propagator, as it is done e.g. in the Brueckner G-matrix derived from a 
realistic NN interaction. The saturation can be obtained if a three-body nucleon is added to the hamiltonian |26l| . 

An alternative technique to determine an effective hamiltonian for a model space calculation is based on a unitary 
transformation of the hamiltonian. It has been developed by Suzuki |27l| and leads to an energy- independent, hcrmitian 
effective interaction. The unitary-model-operator approach (UMOA) has also been used to evaluate the ground-state 
properties of finite nuclei [23, 123, ISfl IsJ . 

In the present study we are going to employ the unitary transformation technique to determine an effective interac- 
tion, which corresponds to the Viow-k discussed above. This effective interaction will then be used in self-consistent 
Green's function (SCGF) calculation of infinite nuclear matter. Various groups have recently developed techniques 
to solve the corresponding equations and determine the energy- and momentum-distribution of the single-particle 
strength in a consistent wav[llli32> 33, 34, 35, 36,i37||. Therefore we can study the correlation effects originating from 
Viow-k inside the model space and compare it to the correlations derived from the bare interaction. Furthermore we 
use the unitary transformation technique to determine an effective interaction which accounts for dispersive effects 
missing in the original Viow-k (see discussion above). 

After this introduction we will present the method for evaluating the effective interaction in section 2 and briefiy 
review the basic features of the SCGF approach in section 3. The results of our investigations are presented in section 
4, which is followed up by the conclusions. 

II. EFFECTIVE INTERACTION 

For the definition and evaluation of an effective interaction to be used in a nuclear structure calculation, which 
is restricted to a subspace of the Hilbert space, the so-called model space, we follow the usual notation and define 
a projection operator P, which projects onto this model space. The operator projecting on the complement of this 
subspace is identified by Q and these operators satisfy the usual relations like P + Q = I, P^ — P, Q^ = Q, and 
PQ = = QP. It is the aim of the Unitary Model Operator Approach (UMOA) to define a unitary transformation 
U in such a way, that the transformed Hamiltonian does not couple the P and Q space, i.e. QU~^HUP = 0. 



For a many-body system the resulting Hamiltonian can be evaluated in a cluster expansion, which leads to many- 
body terms. This is very similar to the folded diagram expansion, which has been discussed above. In UMOA studies 
of finite nuclei terms up to three-body clusters have been evaluated|2al23| indicating a convergence of the expansion 
up to this order. 

In the present study we would like to determine an effective two-body interaction and therefore consider two-body 
systems only. We define the effective interaction as 

Veff^U-^iho + Vi2)U~ho, (1) 

with Vi2 representing the bare NN interaction. The operator /iq denotes the one-body part of the two-body system 
and contains the kinetic energy of the interacting particles. This formulation will lead to an effective interaction 
corresponding to Viow-k- Since, however, we want to determine an effective interaction of two nucleons in the medium 
of nuclear matter, we will also consider the possibility to add a single-particle potential to Hq. Note that in any case 
ho commutes with the projection operators P and Q. 

The operator for the unitary transformation U can be expressed as|38j 

U ={1+L0-UJ^){1+UJLU'' +Uj''uj)-^/^ , (2) 

with an operator uj satisfying u; ~ QluP such that ui'^ = 01^2 — q^ Jj^ ^jjg following we will describe how to determine 
the matrix elements of this operator uj. As a first step we solve the two-body eigenvalue equation 

{ho+vi2)\<^k)=Ek\^k). (3) 

This can be done separately for each partial wave of the two-nucleon problem. Partial waves are identified by total 
angular momentum J, spin S and isospin T. The relative momenta are appropriately discretized such that we can 
reduce the eigenvalue problem to a matrix diagonalization problem. Momenta below the cut-off momentum A define 
the P space and will subsequently be denoted by \p) and \p'). Momenta representing the Q space will be labeled by 
\q) and \q'), while states |i), |j), \k) and \l) refer to basis states of the total P + Q space. 

From the eigenstates |$fc) we determine those Np {Np denoting the dimension of the P space) eigenstatcs \^p), 
which have the largest overlap with the P space and determine 

{q\ij\p') = Y,{q\Q\%){^p\p'), (4) 

p=i 

with \ipp) — P|$p) and {ipp\ denoting the biorthogonal state, satisfying 

^{'^k\p) {p\ipk') and ^{p'\ifk){'Pk\p) ^Spy (5) 

p k 

In the next step we solve the eigenvalue problem in the P space 

^''^IXp) = l4\Xp) , (6) 

and use the results to define 

Wp) = —^\Xp), (7) 

^J■p 

which due to the fact that lu = QluP, can be written as 

{qWp) = -Y.{qHp'){p'\Xp)- (8) 

f^P , 
f pt 

Using Eqs. lO - IJSJ) and the representation of IJ in Eq. lO, the matrix elements of the unitary transformation 
operator U can be written 



(^"\V\p') = (p"|(l+C.t,.)-l/2|p') 



^(I + M^)-i/2(p"|Xp)(Xpb'), (9) 

p=i 



{q\U\p') = {q\uj{\+UJ^Uj)-^'^\p') 
Np 

= Y,{l + ^4r'/^^^,{q\^,){xAp), (10) 

{p'\U\q) = -(pVt(i + ^^t)-i/2|^) 
Np 

= -Y.{l + ^^lVnr>{p'\Xp){^p\q), (11) 

p=i 

{q'\U\q) = (g'|(l + c.^t)-l/2|g) 

= Y.{{l + pl)-"'-l}{q'K){^v\q)+5,,,'- (12) 

p=l 

These matrix elements of U can then be used to determine the matrix elements of the effective interaction Veff 
according to Eq.lQJ. They might also be used to define matrix elements of other effective operators. 

III. SELF-CONSISTENT GREEN'S FUNCTION APPROACH 

One of the key quantities within the Self-consistent Green's Function (SCGF) approach is the retarded single- 
particle (sp) Green's function or sp propagator G{k,uj) (see e.g.|33)- Its imaginary part can be used to determine the 
spectral function 

A{k,uj)^-21mG{k,uj + b]). (13) 

The spectral function provides the information about the energy- and momentum-distribution of the single-particle 
strength, i.e. the probability for adding or removing a particle with momentum k and leaving the residual system at 
an excitation energy related to a;. In the limit of the mean- field or quasi-particle approximation the spectral function 
is represented by a (5-function and takes the simple form 

A{k,u;)^2nS{u;~ek), (14) 

with the quasi-particle energy e^ for a particle with momentum k. The sp Green's function can be obtained from the 
solution of the Dyson equation, which reduces for the system of homogeneous infinite matter to a a simple algebraic 
equation 
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2m ^ ' ' 



G{k,u)^l, (15) 



where I](fc,a') denotes the complex self-energy. The self-energy can be decomposed into a generalized Hartree-Fock 
part plus a dispersive contribution 



m,^)=^--ik)-ir a.''^''^'^-:'^\ (16) 

The next step is to obtain the self energy in terms of the in-medium two-body scattering T matrix. It is possible to 
express Iml](fc,a; + ir]) in terms of the retarded T matrix |llll40l l4l| (for clarity, spin- and isospin quantum number 
are suppressed) 



1 f d'^fc' f^°° dw' 
ImE(fc,c^ + ir;) = ^ J J^ J ^ (kk'|Imr(u; + c.' + iry)|kk') 

x[f{Lu') + b{Lu + u;')]A{k',uj'). (17) 



Here and in the following f{uj) and b{uj) denote the Fermi and Bose distribution functions, respectively. These 
functions depend on the chemical potential ^ and the inverse temperature /3 of the system. The in-medium scattering 
matrix T is to be determined as a solution of the integral equation 

d^qd^c, 

x{qq'\Tin + ir,)\pp') , (18) 



(kk'|r(f} + ir/)|pp') = (kk'll-Ipp') + fflff{kk'\V\qc^)G",,{m'M + ^v) 

J (ZTT) 



where 

J_oo 27r J_^ 2n il - uj - uj' + nj 

stands for the two-particle Green's function of two non-interacting but dressed nucleons. The matrix elements of the 
two-body interaction V represent either the bare NN interaction wi2 or the effective interaction Ve//, in which case 
the integrals are cut at the cut-off parameter A. 

The in-medium scattering equation (|18(l can be reduced to a set of one-dimensional integral equations if the two- 
particle Green's function in H19|) is written as a function of the total and relative momenta of the interacting pair of 
nucleons and the usual angle- average approximation is employed (see e.g. [43 for the accuracy of this approximation). 
This leads to integral equations in the usual partial waves, which can be solved very efficiently if the two-body 
interaction is represented in terms of separable interaction terms of a sufficient rank|33l|. 

Finally, we consider the generalized Hartree-Fock contribution to the self-energy in H16|) . which takes the form 

^"^(k) ^\l ^ (k, k'l V |k, k') n(fc'), (20) 



2 J (27r) 
where nik) is the correlated momentum distribution, which is to be calculated from the spectral function by 

<k)^ f ^ ^fiLu)A{k,Lu). (21) 

J-00 27r 

Also the energy per particle. E/A, can be calculated from the spectral function using Koltun's sum rule 

E 1 f d^fc r+°° duji f e 



A p J (27r)3 y_oo 27r 2 \2m 



Lo]A{k,Lo)f{uj). (22) 



Eas. H13|l - 121|l define the so-called T-matrix approach to the SCGF equations. They form a symmetry conser ving 
approach in the sense of [41|, which means that thcrmodynamical relations like the Hughenholtz-Van Hove theorem |33l 
|43 | are obeyed. 

The Brueckner-Hartree-Fock (BHF) approximation, which is very popular in nuclear physics, can be regarded as 
a simple approximation to this T-matrix approach. In the BHF approximation one reduces the spectral function 
A{k,iLi) to the quasiparticle approximation (|14|l . Furthermore one ignores the hole-hole scattering terms in the 
scattering Ea. (|18|) . which means that one replaces 

{l-f{u)-f{J)) ^ (l-/(c.))(l-/(c.')), (23) 

which is the usual Pauli operator (at finite temperature). This reduces the in- medium scattering equation to the 
Bethe-Goldstone equation. The removal of the hole-hole scattering terms leads to real self-energies I](fc,a;) at energies 
uj below the chemical potential, i.e. for the hole states. 

IV. RESULTS AND DISCUSSION 

In the following we discuss results for symmetric nuclear matter obtained from Self-Consistent Greens Function 
(SCGF) calculations. These calculations are either performed in the complete Hilbert space using the bare CD-Bonn 
m interaction or in the model space, which is defined by a cut-off parameter A = 2 fm~^ in the two-body scattering 
equation, employing the corresponding effective interaction Viow-ki which is derived from the CD-Bonn interaction 
using the techniques described in Sect II. We note that using this unitary model operator technique wc were able to 
reproduce the results of the BHF calculations presented in [23 , which used tabulated matrix elements of |23| , with 
good accuracy. The NN interaction has been restricted to partial waves with total angular momentum J less than 6. 

Results for the calculated energy per nuclcon are displayed in Fig. ^ for various densities, which arc labeled by 
the corresponding Fermi momentum kp. The effective interaction Viow-k accounts for a considerable fraction of the 
short-range NN correlations, which are induced by realistic interactions like the CD-Bonn interactions. Therefore, 
already the Hartree-Fock approximation using this Viow-k yields reasonable results for the energies as can be seen 
from the dotted line of Fig. Q] Hartree-Fock calculations using the bare CD-Bonn interaction yield positive energies 
ranging between 2 MeV per nuclcon and 15 MeV per nucleon for the densities considered in this figure. Note that 
the CD-Bonn interaction should be considered as a soft realistic interaction. Interaction models, which arc based on 
local potentials, like the Argonnc interaction [3, yield more repulsive Hartree-Fock energies Q. 
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FIG. 1: (Color online) Binding energy per nucleon for symmetric nuclear matter as function of the Fermi momentum: Results 
of self-consistent T-matrix calculations for the CD-Bonn potential (dashed line), are compared to results of calculations using 
Viow-k with A = 2fm~^ in the Hartree-Fock approximation (dotted line), the self-consistent second order approximation 
(dashed-dotted line) and for the self-consistent T— matrix approximation (solid line) within the model space. 



The inclusion of correlations within the model space yields a substantial decrease of the energy. The self-consistent 
T-matrix approach provides additional attraction ranging between 6 MeV per nucleon at a density of 0.4 po (with 
Po the empirical saturation density) and 3 MeV per nucleon at 2 po- The fixed cut-off parameter A seems to reduce 
the phase-space available for correlations beyond the mean-field approach at higher densities. Therefore the energy 
calculated in the self-consistent T-matrix approach reduces to the Hartree-Fock result at large densities. 

Fig. ^ also displays the energies resulting from a SCGF calculation within the model space, in which the T-matrix 
has been approximated by the corresponding scattering matrix including only terms up to second order in the NN 
interaction V. The results of such second-order calculations in Viow-k are represented by the dashed-dotted line and 
show a very good agreement with the model-space calculations including the full T-matrix. This confirms the validity 
of approaches, which consider correlation effects within the model-space in a perturbative way. 

All these model space calculations using Viow-k, however, fail to reproduce the results of the SCGF calculations, 
which are obtained in the complete space using the bare NN interaction, which are labeled by CD Bonn T-matrix in 
Fig. Q] In particular, the model space calculations yield to attractive energies at high densities and therefore do not 
exhibit a minimum for the energy as a function of density. This confirms the results of the BHF calculations of |23 . 

It has been argued |23| that this overestimate of the binding energy at high densities is due to the fact that Viow-k 
does not account for the quenching of correlation effects, which is due to the Pauli principle and the dispersive effects 
in the single-particle propagator getting more important with increasing density. Therefore we try to account for the 
dispersive quenching effects by adopting the following two-step procedure. 

In a vein similar to the use of a G-matrix within a self-consistent BHF calculation, as a first step we perform BHF 
calculations using Viow-k- The resulting single-particle spectrum is approximated by an effective mass parameteriza- 
tion. This parameterization of the mean field is employed to define the single-particle operator ho , used in Eq. |^ and 
the following equations of Sect. II (see also jS^)- The resulting effective interaction is used again for a BHF calculation 
within the model space, leading to an update of the mean field parameterization. The procedure is repeated until a 
self-consistent result is obtained. Since the mean field parameterization depends on the density, this method yields 
an effective density-dependent interaction, which in the limit of the density p ^ coincides with Viow-k- Therefore 
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FIG. 2: (Color online) Same as Fig.^but for Viow-k{p) calculated at each density 



we call this effective interaction the density dependent Viow-k or in short Viow-k{p)- Such a procedure amounts to 
summing up certain higher order terms in the full many-body problem. 

In a second step this Viow-k{p) is used in SCGF calculations at the corresponding density. Energies resulting from 
such model space calculations using Viow-k(.p) are presented in Fig. [3 The comparison of the various calculations 
within the model space exhibits the same features as discussed above for the original Viow-k- The correlation within 
the model space provide a substantial reduction of the energy as can be seen from the comparison of the self-consistent 
T-matrix approach with the Hartree-Fock results. The approach treating correlations up to second order in Viow~k{p) 
yields energies which are very close to the complete T-matrix approach. 

The density dependence of the effective interaction Viow-k(p) yields a significant improvement for the comparison 
between the model space calculations and the SCGF calculation using the bare CD-Bonn interaction. Note that the 
energy scale has been adjusted going from Fig.^to Fig. [5] The discrepancy remaining at densities above po niight be 
due to the effects of the Pauli quenching, which are not included in Viow-k{p)- These deviations could also originate 
from the simple parameterization of the dispersive quenching in Viow-k{p)- 

Our investigations also provide the possibility to explore the effects of correlations evaluated within the model space 
using the effective interaction Viow-k- We can furthermore compare these correlation effects with the corresponding 
effects determined by the bare interaction in the unrestricted space. As a first example, we discuss the imaginary 
part of the self-energy calculates at the empirical saturation density po for various nucleon momenta p as displayed 
in Fig. |31 The calculations within the model space reproduce the results of the unrestricted calculations with a good 
accuracy in the energy interval for w ranging between 50 MeV below and 50 MeV above the chemical potential p. 
The remaining differences around the Fermi energy can be attributed to the difference in the effective masses obtained 
using the Viow-k and the bare potential |4J|. The agreement between the T-matrix results around uj ~ p using the 
two potentials is improved if one rescales by the ratio of the effective masses. The imaginary part calculated with 
Viow-ki however, is much smaller than the corresponding result obtained for the bare interaction at energies uj — p 
above 100 MeV. Furthermore the model space calculation do not reproduce the imaginary part for energies below the 
chemical potential at momenta k above 400 MeV/c. 

The imaginary part of the self-energy is a very important ingredient for the evaluation of the spectral function 
A{k,uj) and therefore also for the calculation of the occupation probability n{k) (see Eq. ()21(l ). The small values for 
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FIG. 3: (Color online) Imaginary part of the self-energy as a function of the energy lo for various momenta p as indicated in the 
panels (see Eq. ( 1171 '). The results have been determined for the empirical saturation density po; using Viom-k in the T-matrix 
approximation (solid line), using Viow-k in the second order approximation (dashed-dotted line), and employing CD-Bonn 
interaction in the T- matrix approximation (dotted line) . The dashed line in the first panel denotes the results of the T- matrix 
calculation with the CD-Bonn potential rescaled by the ratio of the effective masses at the Fermi momentum obtained with 
the Viow-k and the bare CD-Bonn potential. 



the imaginary part of the self-energy at high momenta k and negative energies w— /i leads to occupation probabilities at 
these momenta, which are much smaller than the corresponding predictions derived from bare realistic NN interactions 
as can be seen from Fig. 0] This missing strength in the prediction of Viow-k at high momenta is accompanied by 
larger occupation probabilities at low momenta. The self-consistent T-matrix approximation using CD-Bonn yields an 
occupation probability at fc = of 0.897, while the corresponding number using Viow-k is 0.920. At this density, the 
calculation including only terms up to second order in Viow-k yields a rather good approximation to the self-consistent 
T-matrix approximation within the model space. 

As a second example we consider the imaginary part of the self-energy calculated at a lower density p — 0.4 x pQ. 
The results displayed in Fig.|Slrefer to nucleons with momentum k = 0. Also at this density we find that the imaginary 
part evaluated with Viow-k drops to zero at large positive energies much faster than the predictions derived from the 
bare interaction (see upper panel on the left in Fig. |SJ). 

It is worth noting, that at this low density the second order approximation is not such a good approximation to 
the full T-matrix approach as it is for the higher densities. Characteristic differences between the dashed-dotted and 
the solid line show up at energies uj close to the chemical potential. In order to trace the origin of these differences we 
display in Fig.[Slthe contributions of various partial waves of NN interaction channels to this imaginary part. It turns 
out that the differences are largest in the '^S'l — ^ Di and the ^Sq channels. This means that the perturbative approach 
is not very successful in those two channels which tend to form quasi-bound states. In these channels all particle- 
particle hole-hole ladders have to be summed up to obtain the pairing solution. Note, that the pairin g so lutions are 
suppressed at higher densities, if the effects of short-range correlations are properly taken into account |45!, "il?! . 

Furthermore we would like to point out that a different scale is used in the two lower panels of Fig. |5| Taking this 
into account it is evident from this figure that the main contribution to the imaginary part of the self-energy, and that 
means the main contribution to the character of the deviation of the spectral function from the mean-field approach 
originates from the NN interaction in the '^S'l —^ Di channel. 



M 

I" 
I 
O 



C 
O 

■^ 

o 

O 



0.1 



0.01 



0.001 



0.0001 



: 




1 


- ' 


1 


' 1 ' 1 


~ 


- 










V 

low-k 


- 


- 










■-■V, ,2"" order 

low-k 

CDBonn T-matrix 


- 


; 












= 


- 










- 


- 






'^\ 






- 


- 






\A 


L 




- 








S\ 














\ 






— 








^ 


\,\ 


— 


- 










\ N 


- 












\ ^ 




- 












- 












V S 














\ s 




z~ 










\ \ 


—z 


- 










\ \ 


- 


- 










\ \ 


- 


- 










V \ 


- 


- 










\ \ 


- 


- 












- 




1 




1 




, \ 1 , 1^^ 









200 400 600 

Momentum k [MeV/c] 



800 



FIG. 4: (Color online) Momentum distribution n(fe) (see Eq. H'Zll ) calculated for nuclear matter at the empirical saturation 
density pQ. Results of the T-matrix approximation within the model space (solid line) are compared to results of the second 
order approximation (dashed-dotted line) and the T-matrix approximation (dotted line) in the unrestricted space. 



CONCLUSIONS 



During the last few years it has become very popular to perform nuclear structure calculations using effective low- 
momentum NN interactions. These Viow-k interactions are based on a realistic model of the NN interaction. They 
are constructed to be different from zero only within a model space defined by a cut-off A in the relative momenta 
of the interacting nucleons. Within this model space they reproduce the NN data of the underlying bare interaction, 
although the many-body solutions may show differences with different starting NN interactions. 

For this study we performed Self-Consistent Greens Functions (SCGF) calculations of symmetric nuclear matter 
employing Viow-k effective interactions as well as the bare CD Bonn interaction they are based on. Special attention 
was paid to the correlations which can be described within this model space as compared to correlations predicted by 
the underlying interaction within the unrestricted space. 

Using a cut-off A = 2 fm~^ we find that the spectral distribution of the single-particle strength in an energy 
window of plus minus 50 MeV around the Fermi energy is rather well reproduced by the calculation using Viow-k- 
The effective interaction Viow-k is softer than typical realistic NN interactions. Therefore for many observables it is 
sufficient to approximate the full in-medium scattering matrix T by the approximation including terms up to second 
order in Viow-k- This justifies the use of the resummed effective interaction in many-body approximations that do 
not the include ladder-diagram resummation. Special attention must be paid to nuclear systems at smaller densities: 
the possible formation of quasi-bound states may require the non-perturbative treatment of the NN scattering in the 
medium. This also has implications for the use of Viow-k in studies of weakly bound nuclear systems. 

The model space approach cannot reproduce correlation effects, which lead to spectral strength at high energies 
and high momenta. For nuclear matter at the empirical saturation density po the momentum distribution is reliably 
predicted up to a momentum of 400 MeV/c. 

The Viow-k approach overestimates the binding energy per nucleon at high densities. Therefore we introduced a 
density-dependent effective interaction Viow-k{p) which we constructed along the same line as the original Viow-k- 
The new effective interaction accounts for a dispersive correction of the single-particle propagator in the medium. 
This improves the behavior of the effective interaction significantly. For densities above po, however, the binding 
energies calculated with Viow-k(.p) are still too large. This might be improved by determining effective three-nucleon 
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FIG. 5; Imaginary part of the self-energy as a function of the energy uu for nucleons with momentum fc = calculated at the 
density p = 0.4 x po- Results of the T- matrix approach (solid line) and the second order approximation (dashed-dotted line) 
within the model space are compared to results obtained in the unrestricted calculation (dotted line). 



forces explicitly from the underlying bare interaction. 
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